Direct sequencing of bacterial and P1 artificial chromosome-nested deletions for identifying position-specific single-nucleotide polymorphisms

A loxP-transposon retrofitting strategy for generating large nested deletions from one end of the insert DNA in bacterial artificial chromosomes and P1 artificial chromosomes was described recently [Chatterjee, P. K. & Coren, J. S. (1997) Nucleic Acids Res. 25, 2205–2212]. In this report, we combine this procedure with direct sequencing of nested-deletion templates by using primers located in the transposon end to illustrate its value for position-specific single-nucleotide polymorphism (SNP) discovery from chosen regions of large insert clones. A simple ampicillin sensitivity screen was developed to facilitate identification and recovery of deletion clones free of transduced transposon plasmid. This directed approach requires minimal DNA sequencing, and no in vitro subclone library generation; positionally oriented SNPs are a consequence of the method. The procedure is used to discover new SNPs as well as physically map those identified from random subcloned libraries or sequence databases. The deletion templates, positioned SNPs, and markers are also used to orient large insert clones into a contig. The deletion clone can serve as a ready resource for future functional genomic studies because each carries a mammalian cell-specific antibiotic resistance gene from the transposon. Furthermore, the technique should be especially applicable to the analysis of genomes for which a full genome sequence or radiation hybrid cell lines are unavailable.

Whole-exome re-sequencing in a family quartet identifies POP1 mutations as the cause of a novel skeletal dysplasia

Recent advances in DNA sequencing have enabled mapping of genes for monogenic traits in families with small pedigrees and even in unrelated cases. We report the identification of disease-causing mutations in a rare, severe, skeletal dysplasia, studying a family of two healthy unrelated parents and two affected children using whole-exome sequencing. The two affected daughters have clinical and radiographic features suggestive of anauxetic dysplasia (OMIM 607095), a rare form of dwarfism caused by mutations of RMRP. However, mutations of RMRP were excluded in this family by direct sequencing. Our studies identified two novel compound heterozygous loss-of-function mutations in POP1, which encodes a core component of the RNase mitochondrial RNA processing (RNase MRP) complex that directly interacts with the RMRP RNA domains that are affected in anauxetic dysplasia. We demonstrate that these mutations impair the integrity and activity of this complex and that they impair cell proliferation, providing likely molecular and cellular mechanisms by which POP1 mutations cause this severe skeletal dysplasia. © 2011 Glazov et al.

Evolution of genome sequencing techniques

The quality and the speed for genome sequencing has advanced at the same time that technology boundaries are stretched. This advancement has been divided so far in three generations. The first-generation methods enabled sequencing of clonal DNA populations. The second-generation massively increased throughput by parallelizing many reactions while the third-generation methods allow direct sequencing of single DNA molecules. The first techniques to sequence DNA were not developed until the mid-1970s, when two distinct sequencing methods were developed almost simultaneously, one by Alan Maxam and Walter Gilbert, and the other one by Frederick Sanger. The first one is a chemical method to cleave DNA at specific points and the second one uses ddNTPs, which synthesizes a copy from the DNA chain template. Nevertheless, both methods generate fragments of varying lengths that are further electrophoresed. Moreover, it is important to say that until the 1990s, the sequencing of DNA was relatively expensive and it was seen as a long process. Besides, using radiolabeled nucleotides also compounded the problem through safety concerns and prevented the automation. Some advancements within the first generation include the replacement of radioactive labels by fluorescent labeled ddNTPs and cycle sequencing with thermostable DNA polymerase, which allows automation and signal amplification, making the process cheaper, safer and faster. Another method is Pyrosequencing, which is based on the “sequencing by synthesis” principle. It differs from Sanger sequencing, in that it relies on the detection of pyrophosphate release on nucleotide incorporation. By the end of the last millennia, parallelization of this method started the Next Generation Sequencing (NGS) with 454 as the first of many methods that can process multiple samples, calling it the 2º generation sequencing. Here electrophoresis was completely eliminated. One of the methods that is sometimes used is SOLiD, based on sequencing by ligation of fluorescently dye-labeled di-base probes which competes to ligate to the sequencing primer. Specificity of the di-base probe is achieved by interrogating every 1st and 2nd base in each ligation reaction. The widely used Solexa/Illumina method uses modified dNTPs containing so called “reversible terminators” which blocks further polymerization. The terminator also contains a fluorescent label, which can be detected by a camera. Now, the previous step towards the third generation was in charge of Ion Torrent, who developed a technique that is based in a method of “sequencing-by-synthesis”. Its main feature is the detection of hydrogen ions that are released during base incorporation. Likewise, the third generation takes into account nanotechnology advancements for the processing of unique DNA molecules to a real time synthesis sequencing system like PacBio; and finally, the NANOPORE, projected since 1995, also uses Nano-sensors forming channels obtained from bacteria that conducts the sample to a sensor that allows the detection of each nucleotide residue in the DNA strand. The advancements in terms of technology that we have nowadays have been so quick, that it makes wonder: ¿How do we imagine the next generation?

Evidence for u.v. induction of CDKN2 mutations in melanoma cell lines

The CDKN2 gene, encoding the cyclin dependent kinase inhibitor p16, is a tumour suppressor gene involved in melanoma and maps to chromosome band 9p22. Mutations or interstitial deletions of this gene have been found both in the germline of familial melanoma cases and somatically in melanoma cell lines. Previous mutation analyses of melanoma cell lines have indicated a high frequency of C:G to T:A transitions, with all of these mutations occurring at dipyrimidine sites. Including three melanoma cell lines carrying tandem CC to TT mutations, the spectrum of mutations so far reported indicates a possible role for u.v. radiation in the mutagenesis of this gene in some tumours. To further examine this hypothesis we have characterised mutations of the CDKN2 gene in 30 melanoma cell lines. Nineteen lines carried complete or partial homozygous deletions of the gene. Of the remaining cell lines, eight were shown by direct sequencing of PCR products from exon 1 and exon 2 to carry a total of nine different mutations of CDKN2. Two cell lines carried tandem CC to TT mutations and a high rate of C:G to T:A transitions was observed. This study provides further evidence for the role of u.v. light in the genesis of melanoma, with one target being the CDKN2 tumour suppressor gene.

Two novel mutations and a previously unreported intronic polymorphism in the NOTCH3 gene

Cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) is a hereditary disease of small vessel caused by mutations in the NOTCH3 gene (NCBI Gene ID: 4854) located on chromosome 19p13.1. NOTCH3 consists of 33 exons which encode a protein of 2321 amino acids. Exons 3 and 4 were found to be mutation hotspots, containing more than 65% of all CADASIL mutations. We performed direct sequencing on an ABI 3130 Genetic Analyser to screen for mutations and polymorphisms on 300 patients who were clinically suspected to have CADASIL. First, exons 3 and 4 were screened in NOTCH3 and if there were no variations found, then extended CADASIL testing (exons 2, 11, 18 and 19) was offered to patients. Here we report two novel non-synonymous mutations identified in the NOTCH3 gene. The first mutation, located in exon 4 was found in a 49-year-old female and causes an alanine to valine amino acid change at position 202 (605C > T). The second mutation, located in exon 11, was found in a 66-year-old female and causes a cysteine to arginine amino acid change at position 579 (1735T > C). We also report a 46-year-old male with a known polymorphism Thr101Thr (rs3815188) and an unreported polymorphism NM_000435.2:c.679+60G>A observed in intron 4 of the NOTCH3 gene. Although Ala202Ala (rs1043994) is a common polymorphism in the NOTCH3 gene, our reported novel mutation (Ala202Val) causes an amino acid change at the same locus. Our other reported mutation (Cys579Arg) correlates well with other known mutations in NOTCH3, as the majority of the CADASIL-associated mutations in NOTCH3 generally occur in the EGF-like (epidermal growth factor-like) repeat domain, causing a change in the number of cysteine residues. The intronic polymorphism NM_000435.2:c.679+60G>A lies close to the intron–exon boundary and may affect the splicing mechanism in the NOTCH3 gene.

No mutations detected in the INSR gene in a chromosome 19p13 linked migraine pedigree

The aim of this study was to investigate through direct sequencing the insulin receptor (INSR) gene in DNA samples from a migraine affected family previously showing linkage to chromosome 19p13 in an attempt to detect disease associated mutations. Migraine is a common debilitating disorder with a significant genetic component. At present, the number and type of genes involved in the common forms of migraine are not clear. The INSR gene on chromosome 19p13.3-13.2 is a gene of interest since a number of single nucleotide polymorphisms (SNPs) located within the gene have been implicated in migraine with (MA) and without aura (MO). Six DNA samples obtained from non-founding migraine affected members of migraine family 1 (MF1) were used in this study. Genomic DNA was sequenced for the INSR gene in exons 1-22 and the promoter region. In the six migraine family member samples, previously reported SNPs were detected within two exonic DNA coding regions of the INSR gene. These SNPs, in exons 13 and 17, do not alter the normal INSR polypeptide sequence. In addition, intron 7 also revealed a DNA base sequence variation. For the 5' untranslated promoter region of the gene, no mutations or polymorphisms were detected. In conclusion, this study detected no INSR mutations in affected members of a chromosome 19 linked migraine pedigree. Hence, migraine linkage to this chromosomal region may involve other candidate genes.

The reduced expression of the HADH2 protein causes X-linked mental retardation, choreoathetosis, and abnormal behavior

Recently, we defined a new syndromic form of X-linked mental retardation in a 4-generation family with a unique clinical phenotype characterized by mild mental retardation, choreoathetosis, and abnormal behavior (MRXS10). Linkage analysis in this family revealed a candidate region of 13.4 Mb between markers DXS1201 and DXS991 on Xp11; therefore, mutation analysis was performed by direct sequencing in most of the 135 annotated genes located in the region. The gene (HADH2) encoding L-3-hydroxyacyl-CoA dehydrogenase II displayed a sequence alteration (c.574 C-->A; p.R192R) in all patients and carrier females that was absent in unaffected male family members and could not be found in 2,500 control X chromosomes, including in those of 500 healthy males. The silent C-->A substitution is located in exon 5 and was shown by western blot to reduce the amount of HADH2 protein by 60%-70% in the patient. Quantitative in vivo and in vitro expression studies revealed a ratio of splicing transcript amounts different from those normally seen in controls. Apparently, the reduced expression of the wild-type fragment, which results in the decreased protein expression, rather than the increased amount of aberrant splicing fragments of the HADH2 gene, is pathogenic. Our data therefore strongly suggest that reduced expression of the HADH2 protein causes MRXS10, a phenotype different from that caused by 2-methyl-3-hydroxybutyryl-CoA dehydrogenase deficiency, which is a neurodegenerative disorder caused by missense mutations in this multifunctional protein.

Cuff tear arthropathy: Evidence of functional variation in pyrophosphate metabolism genes

We investigated the role of two genes, ANKH and TNAP, in patients with cuff tear arthropathy. These genes encode proteins which regulate the extracellular concentration of inorganic pyrophosphate, fluctuations of which can lead to calcium crystal formation. Variants were detected by direct sequencing of DNA and their frequencies compared with healthy controls. The effect of variants on protein function was further studied by in vitro approaches. Variant genotypes were observed more frequently in the cases when compared with controls in ANKH (45% and 20%) and TNAP (32% and 9%). Variants in ANKH altered inorganic pyrophosphate (PPi) concentrations in transfected human chondrocytes. There was a higher mean serum concentration of TNAP detected in female patients compared with normal ranges. Cuff tear arthropathy is associated with variants in ANKH and TNAP that alter extracellular inorganic pyrophosphate concentrations causing calcium crystal deposition. This supports a theory that genetic variants predispose patients to primary crystal deposition which when combined with a massive rotator cuff tear leads to the development of arthritis.

Novel ANKH amino terminus mutation (Pro5Ser) associated with early-onset calcium pyrophosphate disease with associated phosphaturia

This report describes a 32-year-old woman presenting since childhood with progressive calcium pyrophosphate disease (CPPD), characterized by severe arthropathy and chondrocalcinosis involving multiple peripheral joints and intervertebral disks. Because ANKH mutations have been previously described in familial CPPD, the proband's DNA was assessed at this locus by direct sequencing of promoter and coding regions and revealed 3 sequence variants in ANKH. Sequences of exon 1 revealed a novel isolated nonsynonymous mutation (c.13 C>T), altering amino acid in codon 5 from proline to serine (CCG>TCG). Sequencing of parental DNA revealed an identical mutation in the proband's father but not the mother. Subsequent clinical evaluation demonstrated extensive chondrocalcinosis and degenerative arthropathy in the proband's father. In summary, we report a novel mutation, not previously described, in ANKH exon 1, wherein serine replaces proline, in a case of early-onset severe CPPD associated with metabolic abnormalities, with similar findings in the proband's father.

Pronounced genetic population structure in a potentially vagile fish species (Pristipomoides multidens, Teleostei; Perciformes; Lutjanidae) from the East Indies triangle

The East Indies triangle, bordered by the Phillipines, Malay Peninsula and New Guinea, has a high level of tropical marine species biodiversity. Pristipomoides multidens is a large, long-lived, fecund snapper species that is distributed throughout the East Indies and Indo-Pacific. Samples were analysed from central and eastern Indonesia and northern Australia to test for genetic discontinuities in population structure. Fish (n = 377) were collected from the Indonesian islands of Bali, Sumbawa, Flores, West Timor, Tanimbar and Tual along with 131 fish from two northern Australian locations (Arafura and Timor Seas) from a previous study. Genetic variation in the control region of the mitochondrial genome was assayed using restriction fragment length polymorphism and direct sequencing. Haplotype diversity was high (0.67-0.82), as was intraspecific sequence divergence (range 0-5.8%). FST between pairs of populations ranged from 0 to 0.2753. Genetic subdivision was apparent on a small spatial scale; FST was 0.16 over 191 km (Bali/Sumbawa) and 0.17 over 491 km (Bali/Flores). Constraints to dispersal that contribute to, and maintain, the observed degree of genetic subdivision are experienced presumably by all life history stages of this tropical marine finfish. The constraints may include (1) little or no movement of eggs or larvae, (2) little or no home range or migratory movement of adults and (3) loss of larval cohorts due to transport of larvae away from suitable habitat by prevailing currents

Genetic variability of Tomato spotted wilt virus in Australia and validation of real time RT-PCR for its detection in single and bulked leaf samples

The potential for large-scale use of a sensitive real time reverse transcription polymerase chain reaction (RT-PCR) assay was evaluated for the detection of Tomato spotted wilt virus (TSWV) in single and bulked leaf samples by comparing its sensitivity with that of DAS-ELISA. Using total RNA extracted with RNeasy® or leaf soak methods, real time RT-PCR detected TSWV in all infected samples collected from 16 horticultural crop species (including flowers, herbs and vegetables), two arable crop species, and four weed species by both assays. In samples in which DAS-ELISA had previously detected TSWV, real time RT-PCR was effective at detecting it in leaf tissues of all 22 plant species tested at a wide range of concentrations. Bulk samples required more robust and extensive extraction methods with real time RT-PCR, but it generally detected one infected sample in 1000 uninfected ones. By contrast, ELISA was less sensitive when used to test bulked samples, once detecting up to 1 infected in 800 samples with pepper but never detecting more than 1 infected in 200 samples in tomato and lettuce. It was also less reliable than real time RT-PCR when used to test samples from parts of the leaf where the virus concentration was low. The genetic variability among Australian isolates of TSWV was small. Direct sequencing of a 587 bp region of the nucleoprotein gene (S RNA) of 29 isolates from diverse crops and geographical locations yielded a maximum of only 4.3% nucleotide sequence difference. Phylogenetic analysis revealed no obvious groupings of isolates according to geographic origin or host species. TSWV isolates, that break TSWV resistance genes in tomato or pepper did not differ significantly in the N gene region studied, indicating that a different region of the virus genome is responsible for this trait.

Genetic population structure of red snappers (Lutjanus malabaricus Bloch & Schneider, 1801 and Lutjanus erythropterus Bloch, 1790) in central and eastern Indonesia and northern Australia

The genetic population structure of red snapper Lutjanus malabaricus and Lutjanus erythropterus in eastern Indonesia and northern Australia was investigated by allozyme electrophoresis and sequence variation in the control region of mtDNA. Samples were collected from eight sites in Indonesia and four sites in northern Australia for both species. A total of 13 allozyme loci were scored. More variable loci were observed in L. malabaricus than in L. erythropterus. Sequence variation in the control region (left domain) of the mitochondrial genome was assessed by RFLP and direct sequencing. MtDNA haplotype diversity was high (L. erythropterus, 0.95 and L. malabaricus, 0.97), as was intraspecific sequence divergence, (L. erythropterus, 0.0-12.5% and L. malabaricus, 0.0-9.5%). The pattern of mtDNA haplotype frequencies grouped both species into two broad fisheries stocks with a genetic boundary either between Kupang and Sape (L. malabaricus) or between Kupang and Australian Timor Sea (L. erythropertus). The allozyme analyses revealed similar boundaries for L. erythropterus. Seven allozymes stocks compared to two mtDNA stocks of L. malabaricus including Ambon, which was not sampled with mtDNA, however, were reported. Possible reasons for differences in discrimination between the methods include: i) increased power of multiple allozyme loci over the single mtDNA locus, ii) insufficient gene sampling in the mtDNA control region and iii) relative evolutionary dynamics of nuclear (allozyme loci) and mitochondrial DNA in these taxa. Allozyme and haplotype data did not distinguish separate stocks among the four Australian locations nor the central Indonesian (Bali and Sape locations) for both L. malabaricus and L. erythropterus.

The essential and non-essential genes of Bovine herpesvirus 1

Bovine herpesvirus 1 (BoHV-1) is an economically important pathogen of cattle associated with respiratory and reproductive disease. To further develop BoHV-1 as a vaccine vector, a study was conducted to identify the essential and non-essential genes required for in vitro viability. Randominsertion mutagenesis utilizing a Tn5 transposition system and targeted gene deletion were employed to construct gene disruption and gene deletion libraries, respectively, of an infectious clone of BoHV-1. Transposon insertion position and confirmation of gene deletion were determined by direct sequencing. The essential or non-essential requirement of either transposed or deleted open reading frames (ORFs) was assessed by transfection of respective BoHV-1 DNA into host cells. Of the 73 recognized ORFs encoded by the BoHV-1 genome, 33 were determined to be essential and 36 to be non-essential for virus viability in cell culture; determining the requirement of the two dual copy ORFs was inconclusive. The majority of ORFs were shown to conform to the in vitro requirements of BoHV-1 homologues encoded by human herpesvirus 1 (HHV-1). However, ORFs encoding glycoprotein K (UL53), regulatory, membrane, tegument and capsid proteins (UL54, UL49.5, UL49, UL35, UL20, UL16 and UL7) were shown to differ in requirement when compared to HHV-1-encoded homologues.

Molecular Genetics of Lactase Deficiencies

Congenital lactase deficiency (CLD) (MIM 223000) is a rare autosomal recessive gastrointestinal disorder characterized by watery diarrhea in infants fed with breast milk or other lactose-containing formulas. The CLD locus was previously assigned by linkage and linkage disequilibrium analyses on 2q21 in 19 Finnish families. In this study, the molecular background of this disorder is reported. The CLD locus was refined in 32 CLD patients in 24 families by using microsatellite and single nucleotide polymorphism (SNP) haplotypes. Mutation analyses were performed by direct sequencing. We identified 5 distinct mutations in the lactase (LCT) gene, encoding the enzyme that hydrolyzes lactose in the intestinal lumen. These findings facilitate genetic testing of CLD in clinical practice and enable genetic counseling. The present data also provide the basis for detailed characterization of the molecular pathogenesis of this disorder. Adult-type hypolactasia (MIM 223100) (lactase non-persistence, lactose intolerance) is an autosomal recessive gastrointestinal condition that is a result of a decline in the activity of lactase in the intestinal lumen after weaning. Adult-type hypolactasia is considered to be a normal phenomenon among mammals and symptoms are remarkably milder than experienced in CLD. Recently, a variant C/T-13910 was shown to associate with the adult-type hypolactasia trait, locating 13.9 kb upstream of the LCT gene. In this study, the functional significance of the C/T-13910 variant was determined by studying the LCT mRNA levels in intestinal biopsy samples in children and adults with different genotypes. RT-PCR followed by solid-phase minisequencing was applied to determine the relative expression levels of the LCT alleles using an informative SNP located in exon 1. In children, the C-13910 allele was observed to be downregulated after five years of age in parallel with lactase enzyme activity. The expression of the LCT mRNA in the intestinal mucosa in individuals with the T-13910 A-22018 alleles was 11.5 times higher than that found in individuals with the C-13910, G-22018 alleles. These findings suggest that the C/T-13910 associated with adult-type hypolactasia is associated with the transcriptional regulation of the LCT gene. The presence of the T-13910 A-22018 allele also showed significant elevation lactase activity. Galactose, the hydrolysing product of the milk sugar lactose, has been hypothesized to be poisonous to ovarian epithelial cells. Hence, consumption of dairy products and lactase persistence has been proposed to be a risk factor for ovarian carcinoma. To investigate whether lactase persistence is related to the risk of ovarian carcinoma the C/T-13910 genotype was determined in a cohort of 782 women with ovarian carcinoma 1331 individuals serving as controls. Lactase persistence did not associate significantly with the risk for ovarian carcinoma in the Finnish, in the Polish or in the Swedish populations. The findings do not support the hypothesis that lactase persistence increases the risk for ovarian carcinoma.

Mutation analysis of TGF-β signaling pathway genes among Finnish patients with primary pulmonary hypertension and hereditary hemorrhagic telangiectasia

Primary pulmonary hypertension (PPH), or according to the recent classification idiopathic pulmonary hypertension (IPAH), is a rare, progressive disease of pulmonary vasculature leading to pulmonary hypertension and right heart failure. Most of the patients are sporadic but in about 6% of cases the disease is familial (FPPH). In 2000 two different groups identified the gene predisposing to PPH. This gene, Bone morphogenetic protein receptor type 2 (BMPR2), encodes a subunit of transforming growth factor β (TGF-β) receptor complex. There is a genetic connection between PPH and hereditary hemorrhagic telangiectasia (HHT), a bleeding disorder characterized by local telangiectasias and sometimes with pulmonary hypertension. In HHT, mutations in ALK1 (activin like kinase type 1) and Endoglin, another members of the TGF-β signaling pathway are found. In this study we identified all of the Finnish PPH patients for the years 1986-1999 using the hospital discharge registries of Finnish university hospitals. During this period we found a total of 59 confirmed PPH patients: 55 sporadic and 4 familial representing 3 different families. In 1999 the prevalence of PPH was 5.8 per million and the annual incidence varied between 0.2-1.3 per million. Among 28 PPH patients studied, heterozygous BMPR2 mutations were found in 12% (3/26) of sporadic patients and in 33% of the PPH families (1/3). All the mutations found were different. Large deletions of BMPR2 were excluded by single-stranded chain polymomorphism analysis. As a candidate gene approach we also studied ALK1, Endoglin, Bone Morphogenetic Receptor Type IA (BMPR1A or ALK3), Mothers Against Decapentaplegic Homolog 4 (SMAD4) and Serotonine Transporter Gene (SLC6A4) using single-strand conformational polymorphism (SSCP) analysis and direct sequencing. Among patients and family members studied, we found two mutations in ALK1 in two unrelated samples. We also identified all the HHT patients treated at the Department of Otorhinolaryngology at Helsinki University Central Hospital between the years of 1990-2005 and 8 of the patients were studied for Endoglin and ALK1 mutations using direct sequencing. A total of seven mutations were found and all the mutations were different. The absence of a founder mutation in the Finnish population in both PPH and HHT was somewhat surprising. This suggests that the mutations of BMPR2, ALK1 and Endoglin are quite young and the older mutations have been lost due to repetitive genetic bottlenecks and/or negative selection. Also, other genes than BMPR2 may be involved in the pathogenesis of PPH. No founder mutations were found in PPH or HHT and thus no simple genetic test is available for diagnostics.